Percolation of partially interdependent networks under targeted attack 
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The study of interdependent networks, and in particular the robustness on networks, has attracted considerable 
attention. Recent studies mainly assume that the dependence is fully interdependent. However, targeted attack 
for partially interdependent networks simultaneously has the characteristics of generality in real world. In 
this letter, the comprehensive percolation of generalized framework of partially interdependent networks under 
targeted attack is analyzed. As a = and a = 1, the percolation law is presented. Especially, when a — b — k, 
pi = p2 — p, qa = Is — q, the first and second lines of phase transition coincide with each other. The 
corresponding phase transition diagram and the critical line between the first and the second phase transition 
are found. We find that the tendency of critical line is monotone decreasing with parameter p\. However, for 
different a, the tendency of critical line is monotone increasing with a. In a larger sense, our findings have 
potential application for designing networks with strong robustness and can regulate the robustness of some 
current networks. 

PACS numbers: 89.75.Hc, 64.60.ah, 89.75.Fb 



S3 

cd 

cd 
•*-> 

: 

r/5 ■ 
o ■ 

<Z3 ■ 
>Y 

^3 : 



> 

00 
00 
l> 

in 

od 

o 



I. INTRODUCTION 



Complex networks have been shown to exist in many differ- 
ent areas in the real world and have been intensively studied 
in recent years. However, almost all network research have 
been focused on properties of a single network component that 
does not interact and depend on other networks j l | 4l3ll . Such 
situations rarely, if ever, occur in reality IU4l - [l7ll . In 2010, 
Buldyrev et al. [15] developed a theoretical framework for 
studying the process of cascading failures in one-to-one cor- 
respondent interdependent networks caused by random initial 
failure of nodes. Surprisingly, they found that a broader de- 
gree distribution increased the vulnerability of interdependent 
networks to random failure, which is in contrast to the be- 
havior of a single network. Partially interdependent networks 
were investigated by Parshani et al. Ilcjll . they presented a 
theoretical framework for studying this case in the same year. 
Their findings highlighted that reducing the coupling strength 
could lead to a change from a first to second order percolation 
transition. In 201 1, Huang et al. lfl7ll mainly demonstrated the 
robustness of fully interdependent networks under targeted at- 
tack. The result implied that interdependent networks were 
difficult to defend. However, when real scenarios are con- 
sidered, simultaneous attack on partially interdependent net- 
works intentionally is more general. 

Motivated by the above, we develop a generalized frame- 
work to comprehensively study the percolation of partially in- 
terdependent networks that suffer targeted attack simultane- 
ously. Likewise, percolation law and condition of the first and 
the second phase transitions have been analyzed for partially 
interdependent networks. Furthermore, by applying general- 
ized framework and analyzing the percolation, as a of W a (A;,) 



varying, the corresponding percolation phase transition and 
critical line between the first and the second phase transition 
are firstly presented. Meanwhile, the robustness of two inter- 
dependent networks can be comprehensively studied by using 
the generalized framework. 



II. THE MODEL 

This model consists of two networks A, B with the num- 
ber of nodes Na,Nb, and within each network, the nodes 
are connected with degree distributions Pa(&) and Pg (k), re- 
spectively. Suppose that the average degree of the network A 
is a and the average degree of the network B is b. In addi- 
tion, a fraction qA of network A nodes depends on the nodes 
in network B and a fraction qs of network B nodes depends 
on the nodes in network A. That is, if node Ai of network A 
depends on node Bj of network B and Bj depends on node 
Ak of network A, then k = i. Consequently, when nodes from 
one network fail, the interdependent nodes from the other net- 
work also fail. This invokes an iterative cascade of failures 
in both networks. A value W a (ki) is assigned to each node, 
which presents the probability that a node i with fcj links be- 
comes inactive by targeted-attack. We focus on the family of 
functions 111 811 : 



W a {h) 



kf 



— oo < a < 



(1) 
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When a^O, nodes are attacked intentionally, while for a — 
0, nodes are removed in random. 

We begin by studying the situation where both networks A 
and B are attacked simultaneously with probability W a (ki) 
(Eq. (1)). Initially, 1 — pi and 1 — p2 fraction of nodes are 
intentionally removed from network A and network B respec- 
tively, pa and pb are defined as the fraction of nodes belong- 
ing to the giant components of networks A and B. The re- 
maining fraction of network A nodes after an initial removal 
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of 1 — pi is ip[ — pi, and the remaining fraction of network 
B nodes after an initial removal of 1 — p 2 is <j)' = p 2 . The 
remaining functional part of network A contains a fraction 
ipi = ip'iPA (V4 ) °f network nodes. Accordingly, the remain- 
ing fraction of network B is (f>[ = p2[l — <7s(l — pa(iP[)pi)}, 
and the fraction of nodes in the giant component of network 
B is (j>i — (fiiPE^'i)- Then the sequence, ijj n and <f> n , of giant 
components, and the sequence i// n and <p' n , of the remaining 
fraction of nodes at each stage of the cascading failures, are 
constructed as follows: 

1p[ =pi,lpi = 1p[pA(lp'i), 

00 =P2,(f>'l =P2[l - qs{l -pA(lp[)Pl)],<Pl = (f)'iPB((t>l), 
V^2 =Pl[l - <L4(1 ~ PB((t>l)P2)},1p2 = ^2PA(lp' 2 ), 
<P' 2 =P2[1 - ~PA(ll4i)Pl)]><h = </>2Pb((/>2), 

tfn =Pli 1 - <1A(1 ~ PB(^' n -l)P2)],1pn = V4pa(V4)> 
<P'n =P2[1 - Qb(1 ~ PA(i>'n)Pl)]^n = finPBi&n)- 

(2) 

To determine the state of system (2) at the end of the cas- 
cading process we look at ip' n , cf>' n at the limit of n — > oo. The 
limit must satisfy the equations tp' n — 4i' n+l , <f>' n — 4>'n+i smce 
eventually the clusters stop fragmenting and the fractions of 
randomly removed nodes at step n and n+ 1 are equal. Denot- 
ing ip' n = x, 4>' n — y, we arrive at a system of two symmetric 
equations: 



0.85 



x =pi[l- -Pb{v)P2)] 
V =P2[1- ?b(1 -PA(aj)Pi)] 



(3) 



For equation (1), as a = 0, Fig. 1 show excellent agreement 
between computer simulations of the cascade failures and the 
numerical results obtained by solving system (2). 



m. ANALYTICAL SOLUTION 

Here we discuss the exact analytical results when a = 
and a = 1. As in refs. 01 E3. HHH, we introduce 
the generating function of the degree distribution Gao(0 = 
EfcP/^fc)^, and the generating function of the associated 
branching process, Gai(0 — G' A0 (^)/G' A0 (1). The fraction 
of nodes that belongs to the giant component after the removal 
of 1 — pi nodes is J2J]: 

Pa (pi) = 1-Gao[1-Pi(1- f A )l (4) 
where f A = /a(pi) satisfies a transcendental equation 

f A = G A i[l-pi(l- f A )}- (5) 

When a = 0, Wq = -h, represents the random removal 
of nodes. For the case of two Erdos-Renyi (ER) lfl9l - l2lll 
networks with average degrees a and b, we can easily get 

Pa{x) = 1 — Ja, Pb(d) = 1 — Jb, and system (3) becomes 



X=pi[l- <L4 +P2<L4(1 - Ib)}, 

y = P2[i - m - /a)]- 



(6) 
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FIG. 1: (color online) (a) Simulation results of the giant component 
of the two fully interdependent ER networks after time t cascading 
failures. For each network TV = 100, 000, a — b = 6, pi = 0.8, 
P2 — 0.7, qA — qB = 1. (b) Simulation results of the giant com- 
ponent of the two partially interdependent ER networks after t cas- 
cading failures. For each network JV = 100, 000, a = b — 6, 
qA = 0.6, qB = 0.65, pi — 0.49, P2 = 0.56. All estimates are the 
results of averaging over 10 realizations, ipoo and (poo, the fraction of 
nodes in the giant components of networks A and B separately, after 
a cascade of failures, the simulations results for fully and partially 
interdependent networks are identical with the theoretical values. 



The fraction of nodes in the giant components of networks A 
and B, at the end of the cascading process are given by: 



-0oo =Pi(l - /a)[1 - QA +p 2 <L4(l ~ Jb)], 

0oo = P2(l - - qB +J3l9s(l - /a)]- 



(7) 
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FIG. 2: (color online) For given pi — 0.6 ~ 0.9, with the strong 
coupling, the system undergoes the first phase transition at which 
<j>oa, the fraction of nodes in the giant component of network B, 
abruptly jumps from a finite value to zero. However, with the weak 
coupling, the system undergoes the second phase transition at which 
<j>oo , gradually approaches to zero. 




FIG. 3: (color online) The percolation phase transition for network 
B with a = b — 3, qs = 0.7, and pi = 0.4 ~ 1.0. The critical 
line, as the boundary between the first and second phase transition, is 
founded. The tendency of critical line is monotone decreasing with 
pi . Below the critical line, the system undergoes a first order phase 
transition. Above the critical line, the system undergoes a second 
order transition. 



And /a, Jb can be expressed as: 

f If K 1 r l + fe(Pl 

JA(JB) = 

QB Pi 

f ^ n, 1 r 1 + Qa(P2 - 1) 

JbUa) = — [ 

qA P2 

f B ^l;Vf A J B = l. 



ln/ £ 



bpiP2{fB - 1) 

In !a j 
apiP2(.fA - 1) ' 



(8) 



In fact, for random attack (a = 0) on partially interdependent 
networks simultaneously, the change between the first phase 
transition and the second phase transition can be obtained by 
regulating the coupling strength q A or qs- Fig. 2 shows that 
reducing the coupling strength lead to a change from the first 
to second phase transition. And for different given p\, the crit- 
ical line is also graphically founded ([Fig. 3]). From Fig. 3, 
we can observe that for a = 0, the tendency of critical line 
is monotone decreasing with increasing parameters p\. And 
for pi = 0.6, pi = 0.7, pi = 0.8 and pi = 0.9, the corre- 
sponding percolation phase transitions are shown separately 
in Figure 3. 

Especially when a = b = k, p\ = P2 = p, qA = q_B = q, 
we arrive to the equations: 



Ia = /b = / 



-i)[i-9+p?(i-/)] j0 < / < l, (9) 



and 

V-oo = 0oo = p(l - e- k ^)[l -q + pq{\ - e"**-)]. (10) 

Eq. (9) can be solved graphically as the intersection of a 
straight line y = f and a curve y = e fe p(/- 1 )[ 1 -9+P9( 1 -/)]. 
When p is small enough the curve increases very slowly and 



does not intersect with the straight line except at / = 1 which 
corresponds to the trivial solution. The condition for the first 
order transition (p = p 1 ) is that the derivatives of the equa- 
tions of system (9) with respect to /: 



l = f[kp\l-q) + 2k(p J yq(l-f)] 



(11) 



And solving system (9) for / — > 1 yields the condition for the 
second order transition (p = p 



kp U {l 



q) 



i. 



(12) 



The analysis of Eqs. (11) and (12) show that the first and 
second order transition lines coincide each other: 



P = P 1 ' = P H 



1 



k(l - q) 



(13) 



The critical values of p c , q c for which the phase transition 
changes from the first order to second order are obtained when 
Eqs. (9), (11) and (12) are satisfied simultaneously, we get the 
critical values are only related with average degree as follow- 
ing: 



Pc = 



k + 1- \j2k + 1 



q<- 



V2fc + 1 + 1 
2k ' 



(14) 



The phase transition and critical line are graphically showed 
respectively in Fig. 4. From Fig. 4, the tendency of critical 
line descend with average degree k. 

When a/0, both of the two partially interdependent net- 
works are initially attacked with probability W a (ki) (Eq. (1)). 
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FIG. 4: (color online) The phase transition is showed for k — FIG. 5: (color online) For given pi — 0.6 ~ 0.9, the system un- 
3, 4, 5, 7, 10. The critical line is also graphically presented by Eq. dergoes the first phase transition and the second phase transition for 
(14). Its tendency is monotone decreasing with increasing average k. network B with strong and weak coupling. 



Initially, 1— pi and 1 — p-2 fraction of nodes are intentionally 
removed from network A and network B respectively. The 
generating functions will be defined for network A while simi- 
lar equations describe networks. According to Eq. (1), 1— p\ 
fraction of nodes are removed from network A but before the 
links of the remaining nodes which connect to the removed 
nodes are removed. The generating function of the nodes left 
in network A before removing the links to the removed nodes 

is [nHHm 



Gbo(1 — P2 +P2x) — Gbc{x) and from Eq. (16) , we can get 



Gao(x) = G A b{l + —a;), 

Pi Pi 

Gbo(x) = G b( ,(1-- + -z). 

P2 P2 



(17) 



Simplified forms for GAb{x), Gac{x) and Gao{x) from Eqs. 
(15), (16) and (17) exist when a = 1, 



G Ab (x) = J2 p Pl fa)** 1 = - E p (ki)hfx^ , (15) G ^(x) =^J2 p ( k i) h i^ kl = -GaoM, (18) 

P 1 r. k l 



where the new degree distribution of the remaining frac- 
tion pi of nodes P pi (ki) = ^-P(fci)/iJ 1 , and G a (x) = 

P{k\)x kl , hi = G~ 1 (pi). The generating function of 
the new distribution of nodes left in network A after the links 
to the removed nodes are also removed is 



Gac{x) = G A b{l - Pi +Pix), 



(16) 



is the fraction 



where fi - ElEAiMSl)} - P(fci)fci/ h 1 
wneiepi - - s* 1 P(fci)fc 1 

of the original links that connect to the nodes left, (ki) is the 
average degree of the original network A, (hi (pi)) is the av- 
erage degree of remaining nodes before the links that are dis- 
connected are removed. Then we can find equivalent networks 
A' and B' with generating functions Gao(x) and Gbo(x), 
such that after simultaneous random attack with 1 — p% and 
1 — p2 fractions of nodes, the new generating functions of 
nodes left in A' and B 1 are the same as Ga c (x) and Gbc{x). 
That is, the simultaneous targeted-attack problem on interde- 
pendent networks A and B can be solved as simultaneous 
random-attack problem on interdependent networks A' and 
B' . By solving the equations Gao(1 — Pi + Pix) = Gac{x), 



G Ac {x) = — G A0 {hi{l-pi+Pix)), (19) 
Pi 



G A0 (x) = -G A0 (hi(l - V - + ^x)), 
Pi Pi Pi 



(20) 



where Gao(x) is the original generating function of the net- 
work A, hi = G A l{ Pl ),p! = c^ffi fti- ForERdlH 
networks, we can also get pa(x) = 1 — /a, Pb(v) = 1 — Ib, 
where pa (x) = 1-Gao[1-x{1- Ja)], I a = Gai[1-x(1- 
Ja)], and system (3) becomes 



x = pi[l- qA+P2qA{l - Ib)], 
y =P2[l - qB - /a)]- 



(21) 



The fraction of nodes in the giant components of networks A 
and B, respectively, at the end of the cascade process are then 
given by: 



V>oo =Pl(l - /a)[1 - <L4 +P2<L4(1 - Ib)}, 
<t>oo =P2(1 - /b)[1 - qB +Pl<?B(l - /a)}- 



(22) 
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FIG. 6: (color online) The percolation phase transition for network 
B with a = 3, b = 4, q B = 0.7 and p 1 = 0.6 ~ 1.0. The cor- 
responding phase transition lines are showed. Also, the critical line 
is graphically found, and its tendency is monotone decreasing with 
pi. Below the critical line, the system undergoes a first order phase 
transition. As we approach the critical point, tends to 0. Above 
the critical line, the system undergoes a second order transition. 



And /a, fs can be expressed as: 

f (f s 1 r l + g B (pi -1) 

qs pi 
fA^i;VfA,f B = i, 

1 , l + gA(pa-l) 
q_A P2 



ln/i 



bpip 2 hl{f B - 1) 

In /a j 
ap\Pih\{] A - 1) 



(23) 



/ b 71;V/a,/b = 1. 



Likewise, for targeted attack (a = 1) on partially interdepen- 
dent networks simultaneously, the phase transition changes 
from a first order to a second order percolation transition, as 
the coupling strength q A is reduced (Fig. 5). From Fig. 6, 
the phase transition lines are graphically presented, and the 
critical line descends with increasing parameters p\. 

Especially as a = b = k, p\ — p 2 = p, q A = qs = q, we 
obtain the following equations: 



Ia = /b = / 



,-kph 2 (l-f)[l- q +p q (l-f)] 



and 



)[l-q+pq(l-e 



— kh O 



(24) 



)], (25) 



where h = ^ 



1. The condition for the first order transition 
(p = p 1 ) is that the derivatives of system (24) with respect to 
/: 



1 = f[kp 1 h\l-q) + 2k{p 1 Yqh z {l- /)],0 < / < 1. (26) 

And solving system (24) for / — > 1 yields the condition for 
the second order transition (p = p 11 ): 



FIG. 7: (color online) The corresponding phase transition is showed 
for k = 3, 4, 5, 8, 10. The critical line is also graphically presented 
from simulation. The tendency of critical line is found monotone 
decreasing with average degree k 




FIG. 8: (color online) For given a — 2, 1.5, —2, —3, the system 
undergoes the first phase transition and the second phase transition 
for network A with strong and weak coupling. 



From Eqs. (26) and (27), we have also found that the first 
order transition line coincides with the second order transition 
line (p=p I =p 11 ), the exact condition of phase transition is as 
follows: 



Mp 1V2 1 



(28) 



kp u {\ - q)h z 



(27) 



From Fig. 7, the simulation of phase transition is showed 
for Eq. (28). Meanwhile, the critical line is also by simula- 
tion. From Fig. 7, the tendency of critical line is monotone 
decreasing with average degree k. 

Both networks A and B are attacked simultaneously with 
probability W a (ki), when a > 0, nodes with higher de- 
gree are more vulnerable and those nodes are intentionally 
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FIG. 9: (color online) The percolation phase transition for net- 
work A with a — 3, b — 4, qs = 0.7, p\ — 0.8 and a = 
—3, —2, —1, 0.5, 1.5, 2. The corresponding phase transitions lines 
are showed. Also, the critical line is graphically found, and its ten- 
dency is monotone increasing with a. Below the critical line, the 
system undergoes a first order phase transition. As we approach the 
critical point, tjjoo tends to 0. Above the critical line, the system un- 
dergoes a second order transition. 



attacked, when a < node with higher degree have lower 
probability to fail. For different values of a, Fig. 8 reflect 
the relationship between the phase transition and the coupling 
strength. Fig. 9 show the phase transition lines for different 
given a > respectively. The critical line is also founded and 
its tendency is monotone increasing with a. 



IV. CONCLUSIONS 

For interdependent networks in real scenario, two factors 
are necessary to be considered: partial coupling and targeted 
attack. A general framework is proposed to investigate the 
percolation of partially interdependent networks that suffer 
targeted-attack simultaneously. The percolation of partially 
interdependent networks under targeted attack is comprehen- 
sively analyzed. As a — and a = 1, the percolation law is 
described detailedly . Especially, for a = b = k, p\ = p-2 = p, 
<1a = qs = q, the first and second lines of phase transition 
coincide with each other. And, the tendency of critical line 
monotone decreasing with average degree k. we show both 
analytically and numerically that reducing the coupling be- 
tween the networks leads to a change from a first to a second 
order phase transition at a critical line. The tendency of criti- 
cal line is monotone decreasing with parameter p\ . However, 
for different a, the percolation phase transition is also graph- 
ically demonstrated and critical line is monotone increasing 
with a. Therefore, our finding should be worth considering in 
designing robust network considering. 
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